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The tunneling spectrum of an electron and a hole in a superlattice of NS junctions is computed 
using the BTK approach and the transfer matrix method. It shows sharp resonances at some energies 
above the superconducting gap. The sharper the resonance is the more layers the superlattice 
has. We find for the first time a mechanism to balance the incident and outgoing currents on the 
superlattice by averaging over the phase between the incident electron and the incident hole. This 
mechanism is more natural and physical than those in literatures. 
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I. INTRODUCTION 

The Andreev reflection in normal nietal- 
superconductor(NS) junctions has attracted great 
interests of researchers in the field of superconductivity 
in recent years, 1-4]. The tunneling spectrum of electrons 
on a NS junction is sensitive to the superconducting 
gap of the superconductor (SC), thus providing an 
important technique to measure the gaps and gap 
properties of SCs. When an electron tunnels into a 
SC from a metal there appears a hole reflection to 
conserve the electric charge, vifhich is called the Andreev 
reflection, since the incident electron may form a 
Cooper pair in the SC. A commonly used theoretical 
method to investigate the Andreev reflection is the 
Blonder-Tinkham-Klapwijk(BTK) approachfl|, which 
takes the interface of a NS junction as a 5{x) potential 
barrier. This theory has been widely and successfully 
applied to systems like NS junctions^], ferromagnet-SC 
junctions 0, etc. Wei, Dong and Xing et al studied the 
tunneling spectrum in metal-SC-metal(NSN) junctions 
using the BTK approach^, Q. They found that the 
incident and outgoing currents in both sides do not 
balance each other |7|. Thus they claimed earlier that 
the BTK approach were not suitable for the Andreev 
reflection in NSN junctions and treated them in the 
Landaucr-Buttiker formalism[0|. Later they proposed 
a mechanism to balance the currents by adjusting the 
chemical potential of the SC inside the NSN junctions^ . 
This mechanism seems reasonable since the unbalance 
of currents will change the charge density and thus the 
Fermi surface of the SC. In a similar idea we proposed 
recently another mechanism that the interface of the NS 
junctions is charged by the unbalanced currents[l|- 

Although these mechanisms successfully balanced the 
incident current and the outgoing current, however, they 
mis-considered the phase between the incident electron 



and hole. We found that different phases just balance 
the incident current and the ongoing current. After av- 
eraging over the phases the incident and outgoing cur- 
rents coincide each other through the whole range of bias 
voltage. The charge density, the Fermi surface, and the 
interface charge state won't change at all. This provides 
us a more natural and physical mechanism for the theory 
of multiple NS junctions. In the present paper we study 
the superlattice of multiple NS junctions using this new 
mechanism. 



II. FORMALISM 

A NS junction superlattice is shown in Fig.l, where 
metals and SCs with width L^v and Lg are sandwiched 
periodically but the outmost layers are metals working 
as two leads. A thin insulating interface exists in every 
junction to be treated as a 6{x) potential barrier in the 
BTK approach. An electron and a hole are incident into 
the superlattice under a bias voltage on both sides. There 
is no incident electron in the right metal, thus &„+i = 
and in the left metal no incident hole di = 0. Due to the 
potential barriers the electron and the hole in one layer 
are partially reflected and partially transmitted into the 
next layer. In the SC only a pseudo-particle can propa- 
gate above the energy gap A. When the bias voltage is 
smaller than A only Cooper pairs can do in the SC. This 
causes a hole reflected from the interface, which is called 
the Andreev reflection in literatures. 

The tunneling process of an electron and a hole in the 
NS junction superlattice is described by the Bogoliubov- 
de Gennes (BdG) eg nation 

^.d_(u{v,t)\_(H, A \(u{v,t)\ 

Ho = -^\/' + V{r)~f, (2) 
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where V^(r) = U d{x) on the interfaces. 

When the bias voltage satisfles eV > A an incident 
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FIG. 1: (a) a NS junction superlattice with an incident elec- 
tron and an incident hole under a bias voltage. The metal 
and the SC are colored by blue and yellow, respectively. An 
electron is incident from the left metal and partially tunnels 
into the adjacent SC and is partially reflected. In the right 
side there is no incident electron, 6„+i — and in the left 
side di = 0. (b) The fermi surfaces and energy bands of met- 
als and superconductors. Except for the two metal leads the 
metals and SCs have bias voltage. 



particle in metals tunnels into the SC through the in- 
terfaces and propagate as a quasi-particle inside the SC. 
It is then partially reflected by the following interface 
and partially tunnels into the following metal. The wave 
functions of the electron and the hole inside the zth metal 
layer and the zth SC layer are written as 



N 



Jk+x 



V', 



SC 



-ik-x 

u+ 



(3) 



(4) 



wher 



and 



are the quasi-eletron (e^ > 0) 



v+ J V 

and quasi-hole (e^ < 0) wave functions with the same 
energy E in the SC layers, respectively. They are given 
by 
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The wave vectors are determined by 



E=eV = eV 
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k±a/kF = ± E/fi, 1 < i < n + I 
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When eV < A, q± become complex numbers, so that 
Cooper pairs appear and the wave functions in the SC 
become damping traveling waves. In this case one has 



e, = ±*v/|A|2-£;2 
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(14) 
(15) 



The coefficients a,b,c,d,e, f,g,h in ([3]) to (|3]) are ob- 
tained from the boundary conditions of wave functions 
on the interfaces, where the wave functions are continu- 
ous but due to the d{x) barriers the first derivatives of 
them are not, i.e., 

MXr)^Mx^d (16) 

^b[{X,)-4,'^{X,) + ZkpMXr) = (17) 

where 1 and 2 denote the adjacent metal and SC lay- 
ers, and Z = and Xi is the position of the inter- 
face. These boundary conditions give the following ma- 
trix equations 

(ai,6j,Ci,d^)^ = Mi{Xi){ei,fi,gi,hi)^ (18) 
{eiji,gi,hi)'^ = M2{Yi){ai+i,hi+i,Ci+i,di+i)'^ (19) 

for i — l,2,...,n, where Mi{Xi) and M2{Yi) are called 
the transfer matrices at interface Xi between the metal 
and SC layers and interface Yi between the SC and metal 
layers. Using this matrices we find the relation between 
the coefficients of the incident and outgoing particles in 
the following form 



\{^h{Xi)^h{Y^) 




(20) 



Assigning oi = l,c?i = 0, 6„+i = 0, and c„+i = expiicj)) 
for an incident electron and an incident hole with a phase 
difference 4>, other coefficients can be obtained relative to 
(j) from this matrix equation. These coefficients determine 
the the incident and outgoing current densities and the 
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differential electric conductance [3]. 

Ji ^ [fc+(l- |6i(et/)n + fc_|ci(eT/)p] 
[fc_(l - |d„+i(eF)|2) + fc+|a„+i(et/)|2] 

'^^'^^^^ --l-MeV)f + \c,ieV)f 
|rf„+i(eF)p + |a„+i(e\^)|2 
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In general the incident current and the outgoing cur- 
rent do not match each other in a NS junction superlat- 
tice even in NSN junctions as pointed by Xing et alQ. 
The key idea of this work is that the currents depend on 
the phase difference (j). After averaging over the phase 
the incident current just balances the outgoing current 
and so does the differential conductance. 
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FIG. 2: (a) Probabilities of the reflected electron, the outgo- 
ing electron, the Andreev reflection, and the hole reflection in 
NSN junctions. The total probability of these terms is exactly 
equal to 2. (b) The incident current and the outgoing current 
Ji , Jn+i at different phase 0. Parameters in the computation 
are set to be ^ = 0.5, A = 0.01fj,,Z = l,kFL = 1000, eV = 
2A. 



III. RESULTS OF COMPUTATION 

The phase difference has strong effect on the cur- 
rents on the NS superlattice, as shown in Fig. 2 for a 2- 
NS-junction lattice. It is found that the probabilities of 
the reflected electron |6ip, the outgoing electron |a„+ip, 
the Andreev reflection |cip and the hole reflection |(i„+ip 
vary periodically with the the phase difference (j), and so 
do the incident and outgoing currents. The total prob- 
ability is always conserved at different phase differences 
kn+iP + l^iP + |cip + |rfri+iP = 2. However, the inci- 
dent and outgoing currents do not match each other at 
different phase differences (except two points). This re- 
sult has never been reported in literatures and no one 
noticed that these two currents oscillate about the phase 
difference. Xing et al simply averaged the two currents 
from only in-phase incidences of an electron and a hole 
to balance the two currents |8j- One of the present au- 
thors assumed in an earlier paper a charge accumulation 
on the interfaces to conserve the currents [3]. They all 
neglected the effect of the phase difference on the cur- 
rents. An important result of the present computation is 
that the incident and outgoing currents have almost the 
same average value for different 0. This indicates that a 
great number of electrons and holes with random phase 
differences incident on the two leads of the NS junction 
superlattice will balance the two currents naturally. This 
provides a new but final mechanism for the balance of 
currents on the NS junction superlattice. 

Using the above mechanism the bias voltage is scanned 
for the currents as shown in Fig. 3. It is seen that the cur- 
rents from in-phase incidences of an electron and a hole 
obviously deviate from each other above the supercon- 
ducting gap as reported by earlier wor ksjig]. After av- 
eraging over the phase difference the currents meet each 
other perfectly except small deviations at the three peaks 
of the currents. Therefore, this averaging process success- 
fully recovers the balance of currents in a NS superlat- 
tice. The differential conductance for different interface 
barriers under this phase-averaging mechanism is shown 
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FIG. 3: (a) Incident and outgoing currents in NSN junctions 
with an in-phase incidence of an electron and a hole, (b) The 
currents averaged over the phase 4>. Parameters see Fig.2. 



in Fig. 4(a). For smaller and smaller barrier potentials 
the differential conductance approaches the saturating 
value, i.e. fully conducting NS junctions. As the bar- 
rier increases the differential conductance shows stronger 
and stronger tunneling effect. The first resonance peak 
takes place slightly higher than the superconducting gap. 
This provides a technique for measuring superconduct- 
ing gaps. The other two resonance peaks occur at about 
eV = 1.35A, 1.8A. These resonances need experimental 
verification. For a comparison to Dong's work 8j, where 
only in-phase incidence of electrons and holes is consid- 
ered, Fig. 4(b) shows the differential conductance using 
the phase averaging method and Dong's parameters. For 
weak barriers the present result is quite different from 
Dong's result, but for strong barriers the present result 
approaches to Dong's result. 

Next superlattices of 2,4,8 and 16 NS junctions are 
studied in the above phase averaging mechanism. The 
differential conductances of them with barriers Z = 0.5 
and Z = A are shown in Fig.5(a,b). It is seen that in 
these two cases the NS junction superlattices have some 
sharp zero conductance resonances. Especially those of 
the superlattice of 16 junctions has a full width at half 
maximum of 0.012 A, as seen in the inset of (b), which 
is equivalent to a voltage of about 0.5mV. This effect 
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FIG. 4: Differential conductance for different interface barri- 
ers, (a) parameters see Fig.2; (b) A comparison to Fig. 6 of 
Dong's work[^ with parameters /i — 0.5, A — 0.001/^, kpL = 
5000. 



metal-SC junctions requires a thickness of the SC layer 
above the minimum value. 

In summary, the tunneling spectrum of an electron and 
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FIG. 5: Differential conductances of superlattices of 2,4,8, 
and 16 NS junctions with (a) Z = 0.5 and (b) Z=4.0. Other 
parameters see Fig.2. 

may find applications such as sensitive millivolt electronic 
switches, where a weak voltage increase induces a strong 
current response. 

Finally it should be pointed out that the first reso- 
nance peak of the differential conductance does not oc- 
cur exictly at voltages corresponding to the energy gap 
of the SC at small values of kfL. As shown in Fig. 6 
until kfL — 1200 the first resonance peak moves to the 
gap position. This indicates that the measuring density 
of states using the technique of tunneling spectrum of 



FIG. 6: Differential conductances of multiple NS junction 
lattices with different width kfL of superconducting layers. 
Other parameters see Fig.2. 

a hole on a superlattice of NS junctions is computed us- 
ing the BTK approach. It shows a more abundant struc- 
ture compared to that on a single NS junction, such as 
the sharp resonance above the gap. In particular, the 
sharper the resonance is the more layers the superlattice 
has. This effect may find applications in the electronic 
industry in the future. We find for the first time a mecha- 
nism to balance the incident and outgoing currents on the 
superlattice by averaging over the phase between the inci- 
dent electron and the incident hole. Compared to other 
mechanisms such as the adjusting of the Fermi surface 
and the charge accumulation the present mechanism is 
much more natural and physical. 
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